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(!)  First  Strong  Motion  Records  from  a Central  or  Fast orn 


United  States  Earthquake.  (Fletcher  and  Anderson) 

Two  strong  motion  accelerographs  were  Installed  at 
Blue  Mountain  Lake,  New  York  (BML)  shortly  after  the 
initiation  of  a swarm-type  sequence  of  earthquakes  in 
mid-July  of  1973.  Three  records  rich  in  frequencies  as 
high  as  35  Hz  were  obtained  from  two  of  the  largest  events 
(magnitude  2.7  and  2.6),  the  foci  of  which  are  between 
1.0  and  *1  km  from  the  accelerograph  stations.  The  max- 
imum acceleration  recorded  v/as  +.03*1  g at  ?5  Hz.  These 
are  the  first  strong-motion  records  obtained  from  east  of 
the  Rocky  Mountains.  Since  strong-motion  results  from 
the  western  United  States  may  not  be  applicable  to  the 
eastern  two-thirds  of  the  U.S.,  these  new  results  have 
important  implications  for  the  design  of  critical  facil- 
ities such  as  nuclear  power  plants,  liquified  gas  storage 
tanks,  and  large  dams. 

A paper  on  this  material  has  been  published  (see 
page  4). 

(2)  Source  Properties  of  a Blue  Mountain  Lake  Earthquake. 

(Anderson  and  Fletcher) 

An  accelerogram  obtained  at  Blue  Mt . Lake,  New  York 
is  remarkable  for  the  simplicity  of  its  S-wave  pulse. 

This  results  from  1)  a nearly  complete  absence  of  scatter- 
ing and  reflections  as  second  arrivals  on  the  accelerogram 
and  2)  a very  elementary  earthquake  source.  The  ear';h- 
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quake  identified  with  this  accelerogram  had  a magnitude 
mb  = ?*2  and  a hypocentral  distance  of  about  1 km  from 
the  accelerometer.  Spectral  analysis  of  the  S-wave  indi- 
cates the  earthquake  had  a moment  of  8 x 101B  dyne-cm, 
and  a source  radius  of  20  to  'JO  m.  When  the  accelerogram 
is  integrated  to  obtain  displacement,  there  is  a step 
offset  of  about  5 microns  associated  with  a near  field 
component  of  the  S-wave  pulse.  The  S-wave,  including 
the  step  offset,  can  be  matched  in  remarkable  detail  by 
a dislocation  model  with  a moment  of  8.'l  x 101B  dyne-cm. 

A paper  on  this  material  has  been  published  (see 
page  4). 

( 3 ^ Comparison  of  Strong  Ground  Motion  from  Several 

PA?locat lop  Models.  (Anderson  and  Richards) 

To  examine  the  effects  of  different  earthquake  sources 
ground  motion  near  the  rupture  for  several  dislocation 
models  has  been  calculated.  Quite  different  dislocation 
models  of  the  r \pture  can  give  very  similar  displacements 
at  stations  located  only  one  fault  dimension  from  the  fault 
in  an  infinite  homogeneous  elastic  space.  In  particular: 

1)  models  with  fault  motions  and  rupture  geometry  quite 
different  from  Haskell’s  propagating  ramp  can  often  have 
near-field  motions  very  similar  to  those  from  a ramp  model; 
and  2)  quite  similar  near-field  motions  are  found,  for 
different  ramp  models  in  which  the  rupture  velocity  and 


rise  time  are  varied  together  over  rather  large  ranges 
(e.g.,  a factor  of  three  in  rise  time). 

V/e  infer  that  there  is  considerable  ambiguity  in 
interpreting  near-field  displacement  data  In  terms  of  a 
model  of  fault  motion.  The  ambiguity  may  be  reduced  by 
using  acceleration  data,  although  tills  entails  a consid- 
erable increase  in  computer  time  for  solving  forward 
problems . 

A paper  on  tills  material  has  been  published  (see 
page  H ) . 

( ) The  Dynamic  Field  of  a Growing  Plane  Elliptical 

Shear  Crack.  (Richards) 

The  radiation  for  a three-dimensional  problem  of 
brittle  fracture  is  investigated.  A crack  is  presumed 
to  nucleate  at  a point  in  an  infinite  pre-stressed 
elastic  medium,  and  the  crack  subsequently  grows  steadily 
with  subsonic  rupture  velocities,  maintaining  the  shape 
of  an  ellipse.  Shear  stresses  are  relieved  by  the  crack, 
and  exact  solutions  are  derived  for  the  acceleration  and 
stress-rate  (at  every  point  of  the  medium)  in  terms  of 
single  Integrals  and  algebraic  expressions.  The  solutions 
are  evaluated  analytically  at  wavefronts  and  singularities, 
and  numerically,  at  different  points  in  the  medium,  for 
different  growth  rates  of  the  crack. 

A paper  on  this  material  has  been  published  (see 
page  *0- 
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( 5 ) Motions  Near  a Shallow  Rupturing  Fault:  Evaluation 


of  Effects  due  to  the  Free  Surface.  (Anderson) 

A full-size  report  on  this  topic  is  attached  herewith 
as  an  appendix.  Report  consists  of  22  text  pages,  5 
appendix  pages,  2 reference  pages,  17  figure  caption  and 
figure  pages,  and  1 table. 

This  material  has  been  accepted  for  publication  (see 
page  4). 
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MOTIONS  NEAR  A SHALLOW  RUPTURING  FAULT: 

EVALUATION  OF  EFFECTS  DUE  TO  THE  FREE  SURFACE* 

John  G.  Anderson 

Lamont -Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  New  York  10964 
and 

Department  of  Geological  Sciences 
Columbia  University 

ABSTRACT 

The  displacement  of  the  surface  of  a half  space  near  a 
shallow  rupturing  fault  is,  generally,  approximated  poorly 
by  the  method  of  doubling  the  amplitude  calculated  for  the 
same  source  in  an  infinite  space.  To  obtain  this  result, 
the  motions  of  a half  space  were  calculated  using  a Green's 
function  which  is  a solution  to  Lamb's  problem,  and  the 
motions  of  an  infinite  space  were  calculated  using  the 
formulae  of  Haskell.  The  infinite  space  method  does  work 
well  for  two  cases:  1.  A fault  from  which  the  angles  of 

incidence  at  the  station  are  less  than  30  degrees  and  2. 
Horizontal  components  of  motion  dominated  by  SH  waves  (as 

from  a strike  slip  fault)  with  angles  of  incidence  less  than 
60  degrees. 


•Lamont-Doherty  Geological  Observatory  Contribution  No.  0000. 


2. 


INTRODUCTION 

The  surface  notions  resulting  from  a point  source  in 
a naif  spoce  include  four  major  effects  which  arc  absent 
in  an  infinite  space:  The  amplification  of  all  waves;  the 

phase  snilt  of  dV-wavos  incident  at  angles  greater  than 
critical;  the  dP-phase,  which  travels  to  the  surface  as 
an  07-wave  at  the  critical  anrlo  and  is  then  refracted 
horizontally  as  a ?-.:r,ve;  and  Ray  lei,*;’,  waves . These 
ejects  have  all  been  studies  before  from  a point  source 
(e.f*.  hnopoff  ci  ul.  t 1997;  Pekeris  and  hi f sen,  1?97; 
nawasaki  et  al.3  1073),  the  present  objective  is  to  oval- 
uatc  t.ie  ex  feet  of  including*  the  free  surface  when  the 
source  is  a spatially  extended  rupture.  This  allows  a 
critical  evaluation  of  the  assumption,  used  frequently  In 
dislocation  modelim;;,  that  tne  free  surface  may  be  reasonably 
accounted  for  by  douolin^*  tne  amplitude  of  notions  in  an 
infinite  space  (e.£*.  r.anarriori,  1972;  Trifunac,  197*1;  Trifunac 
and  Udwadia,  197*0  Anderson,  197*0. 


r 


;rj?i:oD 

The  displacements  of  the  surface  of  a half  space  are 
calculated  by  applying  the  representation  theorem  (furri-J  n 
and  hnopoff,  lQCdl): 


U±(x,t) 


/ dt’  ffv n 


k(OCuj(c,t»)] 


Here  U±  Is  the  ifc/t  component  of  displacement  of  the  sur- 
face, y the  shear  modulus,  nk  the  kt/z  component  of  the 
normal  to  the  fault,  and  [u^]  the  amplitude  of  the  dis- 
location on  the  fault.  The  integration  is  over  time  and 
the  fault  surfaces.  The  Green’s  function  is 


, where  ic  the  displacement  in 


the  i direction  at  location  x,  time  t due  to  an  impulse 
force  in  the  j direction  at  location  £,  time  t’.  Thus, 
gij,k'  is  the  displacement  due  to  a point  couple,  and 


the  sum  (g 


ij  ' 


is  the  displacement  due  to  a 


double  couple  force,  or,  equivalently , 
dislocation.  It  is  of  the  form 


di_ 

dt2 


h 


ij,k’ 


due  to  a point 


glJ,k' 


4. 


where  h 


ij.lt' 


is  evaluated  by  an  integral  in  the  conpl.  :• 


plane  given  by  Johnson  (197*0.  mhe  differentiation  met 

be  done  numerically.  As  a function  of  tine  t = t-i.  ’ 

o 

h is  a Green's  function  for  a force  which  is  ;:cro  fo: 

1 J , K 

t <0  and  has  magnitude  t for  t >0. 

o oo 

This  series  of  integrations  and  differentiations  was 
handled  as  follows: 

Step  1.  The  Green's  functions  are  computed  (us in.*  a ".O’ 
integration  subroutine)  and  stored  for  an  ;-”rey 
of  points  on  the  fault.  The  differentiation  •••*  f- 


respect  to  tine  is  deferred  until  the  second  stop. 
These  are  computed  only  once  for  each  fault-sin tion 
geometry. 

Step  2.  The  source  tine  function,  for  each  grid  point  Is 

convolved  with  the  corresponding  Green's  function 
to  perform,  the  time  integral.  The  results  from  all 
the  points  on  the  fault  plane  are  weighted  and 
summed  to  perform  the  spatial  integration.  The 
weight  of  a grid  point  is  the  area  of  the  fault  in- 
side a rectangle,  centered  on  the  grid  point,  with 
sides  equal  to  the  grid  spacing.  Thus,  the  sum  of 
the  weights  is  equal  to  the  fault  area.  The  sun 
of  the  convolved  functions  is  differentiated  numeri- 
cally to  obtain  the  desired  displacement  at  the 


station . 


Because  the  spatial  Units  of  integration  are 


a function  of  tine,  it  ic  necessary  to  Justify 
taking  the  tine  differentiation  outside  the 
spatial  integral.  This  is  allowable  because  the 
integration  is  approximated  numerically  by  a sum, 
v/here  the  weight  of  the  Green’s  function  for  each 
grid  point  is  independent  of  tine.  The  errors  due 
to  this  approximation  are  discussed  below.  lo- 
calise the  weights  are  constants,  it  does  not 
matter  whether  the  sum  or  the  numerical  differ- 
entiation is  dene  first. 


In  finding  a dislocation  model  for  an  earthquake,  the 
fault-station  geometry  is  known,  and  displacements  at  the 
station  are  typically  calculated  for  several  trial  models 
of  rupture.  The  two  step  method  is  excellent  in  this 
situation,  because  storing  the  Green’s  functions  as  in 
step  1 is  more  economical  than  recomputing  them. 

The  displacements  from  a dislocation  in  an  irkin’ 
space  are  calculated  as  described  in  indorsor,  and  Richards 
(1975).  The  coordinate  system  for  all  models  is  showr  ,1.r 
Figure  1.  The  F-wave  velocity  n = f.°  !;m/~r»c  and  f-mve 
velocity  f = km/soc  war.  usod  throughout. 

A propagating  ramp  dislocation  tine  function  is  used 
here,  but  any  other  time  function  could  have  been  used. 


In  this  model, 


where  L , L are  the  fault  dimensions,  v is  the  rupture 

velocity,  x is  the  rise  time,  and  D is  the  final  dis- 

m 

placement  cn  the  fault. 
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obtain  the  Green’s 

To  estimate  f . 

s * 


function  is  considered, 
consider  a source  with  an  impulse 


time  function  traveling  with  velocity  v along  a short 


line  of  length  L.  The  radiation  at  a distant  station  is 
a square  wave,  with  duration  T#  and  amplitude  1,  say. 

The  Fourier  amplitude  spectrum  is 


7. 


F (u)  = T gAf?  X 

<2  0 A 

cjT 

where  X = . A model  for  this  source  by  summing  the 

2 

radiation  from  two  stationary  points  at  the  ends  gives  an 
approximate  waveform 

fft(t)  = T0  6 /t  + M + ^o_  6 /t  - 

~ y 2 J 2 l 2 

The  Fourier  amplitude  spectrum  for  this  approximation  is 

F (uj)  = T cos  X. 

3.  0 

At  zero  frequency  both  the  approximate  and  the  exact 
spectra  have  amplitude  T and  zero  slope,  but  as  |x|  in- 
creases they  diverge.  We  introduce  X as  the  upper  limit 

O 

on  values  of  |x|  for  which  F is  a satisfactory  approxi- 
mation  to  F . X should  almost  certainly  be  given  a value 

G S 

less  than  1,  for  at  X = 1,  F„  is  only  about  two-thirds  of 

cl 

Fg.  For  a propagating  fault,  Ben-Menahem  (1961)  showed 

X = uL  (c-  cos  0) 

2c  v 

where  c is  the  wave  propagation  speed  and  6 is  the  angle 
between  the  direction  of  rupture  propagation  and  the  dir- 
ection to  the  receiver.  If  |x|  < X , then  the  frequencies 

f which  are  adequately  represented  are  given  by 


A 


a 
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V 

w < c ' s = f . 

2 7T  tFU  (C.--CGS  0 ) 
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Thus,  a more  dense  grid  spacing  is  needed  to  obtain  dis- 
placements at  backward  angles  fron  the  direction  of  ruptur 
propagation.  For  the  cases  of  this  paper,  we  use  the  sim- 
plified estimate  f = 3 , which  can  be  derived  using 

Trn 

less  than  0.6  for  our  source  receiver  geometries.  For 
other  cases,  particularly  -..'here  0 > n or  for  small  rupture 


velocities,  the  exact  formula  should  oe  used. 

Computer  time  is  in  general  proportional  to  f 3 . ’ o 

powers  of  fg  arise  because  the  di.3tfu.cu  L is  used  as  ‘ ; 


Greer.1 

's  functions 

er.  the 

U1 1 in  c t;| , 

coor  d*  •.r,t 

The 

' hire;  power 

arises 

because  ror  a 

lr."'*or 

fs  the  Green's  function  must  be  calculated  at  a greater 
number  of  points  in  time.  Here,  waves  at  frequency  f„  are 
sampled  at  a rate  of  six  times  per  cycle. 

The  effect  of  fg  may  be  illustrated  by  replacing  the 
Green's  function  in  stop  1 with  the  nonphysical  function 
h = 1 II(t  - R/cO,  where  !i  is  the  Heaviside  step  function 


and  R is  the  distance  between  the  grid  point  on  the  fault 
and  the  receiver.  The  discontinuity  in  h is  typical  of 
realistic  Green's  functions.  A model  in  shown  in  rimurc  P, 
using  this  function  h,  for  three  values  or  f 


mhoso  value's 


9. 


required,  respectively,  7,  180,  and  2303  total  evalutions  of 

h in  space  and  time.  Thus,  each  successively  greater  value 

of  fs  caused  over  »*m  order  of  magnitude  increase  in  computer 

time.  The  lines  beneath  the  time  axis  in  Figure  2 show  the 

time  interval  T = f”1  for  each  value  of  f „ . One  obvious 
s s s 

inaccuracy  is  that  the  computed  waves  do  not  begin  at  the 
theoretical  arrival  time  shown  by  the  arrow.  This  discrepancy 
is  inevitable  for  the  case  with  only  one  grid  point,  ns  that 
point  is  not  at  the  origin  of  rupture.  For  the  other  two 
cases,  where  a grid  point  is  at  the  origin  of  rupture,  the 
discrepancy  originates  in  the  numerical  second  derivative. 

In  all  cases,  the  difference  between  the  time  when  the  computed 
wave  becomes  non-zero  and  the  theoretical  arrival  time  is 
considerably  shorter  than  Ts.  The  waveform  would  not  be 
much  different  for  values  of  f5  larger  than  the  largest 
shown,  where  the  duration  of  the  wave  is  about  four  times  T s . 

The  integration  to  evaluate  the  hrcca's  function  in 
step  1 is  the  final  major  source  of  errors  in  the  no:  rutod 
waveform.  This  is  done  in  a hcmber-  integration  subroutine 

(Uilf,  1967),  which  forms  a series  of  estimates  to  the  inte- 
gral, doubling  the  number  of  points  in  the  integrand  for 

each  successive  estimate.  It  uses  this  series  to  predict 
the  value  of  the  integral,  and  returns  an  answer  when  two 
successive  predictions  differ  by  less  than  a specified  rela- 
tive error  (r,  say).  In  computin'-  the  hrorn’s  '’unction  °or 
a series  of  times  in  step  1,  th'"  answer  e-’ch  M.r''  nould 
have  a random  error,  or  the  answer"  »***.»»  pr»  a ”r t p":.a t i e ^ ] 1 


CO 
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The  half  space  Green's  functions  were  subjected  to  sev 
eral  tests,  to  be  sure  they  are  computed  correctly.  These 
will  now  be  mentioned  briefly  before  preceding  to  some  com- 
parisons of  half  space  and  infinite  space  displacements. 

1*  Displacements  agreed  with  Figures  7,  °,  and  ; of 
Johnson  ( 197*1 ). 
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First  motions  were  in  the  proper  directions. 
Rayleigh  waves  have  retrograde  elliptical  particle 


motion 
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Symmetry  properties:  All  components  of  radiation  from 

several  sources  showed  the  expected  behavior  of  being 
either  even  or  odd,  depending  on  the  source.  In  both 
x and  y. 

5.  The  static  offsets  vrere  in  accord  with  results  of  Sato 
and  Matsuura  (197*0. 

6.  Wave  equation:  Using  reciprocity,  a Green’s  function 

Sfj  k»  will  obey  the  wave  equation  for  a receiver 
fixed  on  the  free  surface  v; hen  the  source  location  is 
varied  to  compute  the  spatial  derivative.  As  pro- 
grammed, the  terms  g in  eq.  (1)  are 

summed  algebraically.  For  P-waves,  g^.  k,  = g^ 

and  computed  displacements  were  shown  to  obey  the  wave 
equation.  The  3-component  lacks  this  symmetry,  and 
could  not  be  subjected  to  this  test. 

The  comparison  in  the  next  section  of  half  space  and 
infinite  space  displacements  from  a point  dislocation  source 
may  also  be  regarded  as  a mutual  check  of  the  half  space 
and  infinite  space  calculations. 

NUMERICAL  EXAMPLES:  POINT  DISLOCATION  SOURCES 

Dynamic  displacements  at  a distance  r from  a point 
dislocation  source  in  an  infinite  space  consist  of  five 
components:  two  far  field  terms  which  decrease  as  r“l, 

two  Intermediate  field  terms  which  decrease  as  r-2 , and 
the  near  field  term  which  decreases  as  r“4  (Haskell,  1969). 
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A comparison  of  motions  from  a point  dislocation  in 
a half  space  and  in  an  infinite  space  serves  two  purposes. 
First,  it  shows  how  well  the  far  field  components  from 
a point  source  in  an  infinite  space  (which  can  be  quickly 
computed)  can  be  corrected  using  plane  wave  theory  to 
obtain  the  corresponding  components  in  the  half  space, 
decor.d,  it  shows  how  the  near  and  intermediate  field  dis- 
placements are  affected  by  the  free  surface. 


Displacements  were  calculated  at  eight  stations, 
chosen  such  that  in  the  half  space  direct  rays  have  anal  ■'s 
of  incidence  (measured  fro;:  the  normal  to  the  surface) 
which  varied  from  10  to  20  degrees  ir.  10  degree  intervals. 
Wc  used  two  sources,  both  on  a vertical  fault:  one  was 

strike  slip  and  the  other  was  dip  slip,  flic  stations 

are  on  the  x axis,  so  that  far  field  radiation  is  cu- 

2 

tirely  S1I  from,  the  strike  slip  source,  and  entirely 
P-3V  from  the  dip-slip  source,  h'ote  that  the  displace- 
ments for  the  Sh  case  cannot  be  calculated  by  the 
source-image  method  because  the  dislocation  source 
causes  P-  and  SV-far  field  motions  elsewhere  in  the 
* 0 plane. 

From  these  calculations,  the  free  surface  amnliricatio 
.was  derived  for  a far  field  body  wave  by  dividing  the  com- 
puted amplitude  for  the  half  space  by  the  corresponding 
amplitude  for  the  infinite  space  case.  Dor  P-v.avcs,  the 
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computed  ratios  arc  shown  in  Figure  , and  for  SV-wavec, 
these  ratios  arc  shown  in  Figure  Z. 

The  half  space  SH  waves  had  the  same  wave  shape  as  the 
whole  space  waves,  and  twice  the  amplitude,  to  within  32. 
This  agrees  well  with  theory,  which  predicts  that  exactly 
twice  the  amplitude  should  be  expected  at  all  angles  of 
incidence . 

For  P and  SV  waves,  the  theoretical  amplification  of 
waves  by  the  free  surface  is  given  by  equations  in  the 
Appendix.  Corresponding  theoretical  curves  are  plotted 
in  Figure  4 for  P-waves  and  Figure  b for  SV-waves,  to- 
gether with  the  ratios  of  the  components  of  motions  de- 
rived from  the  model  calculations.  For  P-waves,  as 
discussed  by  Kawasaki  et  at.  (1373),  the  free  surface  aoes 
not  greatly  affect  the  wave  shape,  and  the  calculations 
agree  well  with  the  theory. 

Because  the  SV  displacements  are  strongly  modified  by 
the  free  surface,  they  are  shown  for  both  the  half  space 
and  the  whole  space  in  Figure  6 for  stations  at  30  to 
80  degrees.  The  station  at  30  degrees  is  not  beyond  the 
critical  angle  (3^.7°),  there  is  no  phase  shift,  and  the 
wave  shapes  are  similar.  For  the  other  stations,  we  must 
actually  apply  the  theoretical  phase  shift  to  determine 
how  well  the  calculations  agree  with  theory.  To  evaluate 
the  effect  of  this  phase  shift,  consider  an  incident  wave- 
form f(t)  with  Fourier  transform  F(u>).  The  Fourier  trans- 
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form  of  the  surface  motion  Is,  using  formula  A-5  or  A-6, 

U(w)  = R e1<SsCn(to)  p(w)>  where  R is  |RXS|  or  |rzs|  , and 

6sgn(w)  is  the  corresponding  phase  shift.  Applying  the 
inverse  transform  to  obtain  the  surface  motion  gives: 


00 

u(t)  =_1  f Rei6sr,n(w)  P(w)elfr>td 
2 ir 


B R -os  6 f(t)  + R sin  5 A(f(t)) 
where  // ( f ( t ) ) = 


00  00 

s I f d“  / f ( t ) sin  u (t  - t)  dr. 


The  function  //(f(t))  is  the  Hilbert  transform  of  f(t), 
discussed  in  more  detail  by  Choy  and  Richards  (1975). 

The  phase  shifted,  amplitude  corrected,  -d  ole 
SV  waves  are  compared  ,-ith  the  half  space  vc.s  1r  -l;«c 
7.  Except  at  4o°  an-  3"°  the  exact  :;avr  shapes  ar*'li- 
tuder.  agree  v:cll  with  these  found  by  applying  the  plane  pave 
correction.  At  40°  the  displacements  shov:  rractioa.il-'  no 
resemblence  to  each  other.  At  50°  there  are  lesser  diff- 
erences, and  another  phase  arrives  before  the  SV-wave . 

This  is  the  SP-phase.  At  tO°,  this  phase  Is  superimposed 
on  the  S-wavc , and  is  perhaps  the  sole  cause  of  the  dis- 
agreement. Although  for  angles  greater  than  50°  the  SP- 
phase  has  moved  out  of  the  time  window  in  Figure  7,  it 


appears  out  to  80°,  with  a shape  similar  to  the  SP-phase 
at  50° , but  decreasing  amplitudes  relative  to  the  S-wave. 

As  at  ^0°  and  50°,  amplitudes  of  the  SP-phase  may  be  larger 
than  the  S-wave,  but  in  Figure  7 it  appears  to  be  depleted 
in  high  frequencies  relative  to  the  S-wave. 

The  second  objective  of  comparing  displacements  from 
a point  dislocation  in  a half  and  an  infinite  space  was  to 
study  the  amplification  of  near  and  intermediate  field 
components.  These  displacements  cannot  be  separated  for 
individual  study  as  were  the  far  field  terms.  Dut  at  times 
when  far  field  terms  are  absent,  the  ratios  of  the  half 
space  to  the  infinite  space  displacements  reveal  the  net 
effect  of  the  free  surface  upon  any  near  field  or  inter- 
mediate field  terms  which  are  present.  This  net  er'‘oe4-  is 

complex,  and  in  general  is  not  well  described  by  the 
approximation  that  the  free  surface  causes  the  amplitudes 

to  double. 

Because  the  stations  are  on  the  x axis  (Figure  1, 

2 

Figure  6 caption),  the  SH  case  causes  only  U displacements 
and  the  P-SV  case  causes  only  U (horizontal)  and  U 

2 3 

(vertical)  displacements.  In  speaking  of  the  near  and 
intermediate  field  terms.  It  is  better  to  refer  only  to 
the  U^,  and  components  to  remove  the  connotation  of 
phases  travelling  only  with  the  P-  and  S-velocities . 

Immediately  after  the  P-wave,  the  amplitude  ratio 
for  each  component  resembles  the  amplitude  ratio  for  the 
horizontal  component  of  the  P-wave  (Figure  *0 . Just  before 
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the  S-wave,  the  ratio  lor  each  component,  where  It  could 
be  measured,  had  Increased  to  a value  generally  In  the 

ranre  of  2.5  to  2.8.  This  ratio  could  not  !-o  rcisur-'’ 
the  H2  and  V components  Incident  at  angles  -r^atm  '’'an 
decrees  bececse  In  those  cases  the  rp-voyo  n-d  ~v  „ „ 
shift  of  the  P— wave  also  cause  •»  ♦ >,ri  > nio 

space  before  the  then***tical  p_--nV^  . V!o  h,t 

The  static  offset  ratios  for  all  three  components,  shown 
In  Figure  8,  generally  differ  from  2.0.  These  ratios,  like 
the  ratios  of  far  field  components,  are  Independent  of  the 
distance  between  the  source  and  the  station.  For  angles  of 
Incidence  less  than  30  degrees  these  ratios  differ  consider- 
ably from  the  amplification  ratios  of  far  field  body  waves. 
Thus,  It  appears  that  even  when  the  dynamic  displacements 

of  the  surface  of  a half  space  can  be  approximated  by  twice 
the  Infinite  space  displacements,  the  static  displacements 
cannot  be  reliably  approximated  In  this  way. 

In  summary,  the  amplification  and  phase  shift  of  the 
far  field  components  of  displacement  close  to  a point  dis- 
location can  be  understood  well  by  applying  plane  wave 
theory.  The  theoretical  amplification  of  plane  waves  does 
not  apply  to  the  static  offset  or  to  the  amplification  of 
near  field  components  of  dynamic  displacements,  and  it  can- 
not, of  course,  explain  the  SP-phase  or  the  Rayleigh  wave. 
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NUMERICAL 


EXAMPLES : 


EXTENDED  DISLOCATION  SOURCES 


For  a snail  earthquake  (resembling  a point  source), 
Figures  4-7  show  that  doubling  infinite  space  notions  is 
inadequate  for  stations  with  an  angle  of  incidence  of  over 
30  degrees  if  there  is  any  SV  notion,  and  nearly  always 
inadequate  if  the  angle  of  incidence  is  greater  than  79 


degrees.  For  an  extended  source,  however,  the  contribution 
from  each  Green's  function  is  only  a snail  part,  of  the 
total  notion  at  the  station.  Therefore  Figures  9- 13  were 
drawn  to  study  half  space  and  -whole  space  notions  for 
extended  vertical  faults.  The  fault  motion  is  dip  slip 
for  the  models  in  Figures  9 and  19,  so  that  SV  notion  dom- 
inates, ana  strike  slip  for  the  models  in  Figures  11-13,  so 
that  SH  motion  dominates.  The  angles  of  incidence  from  the 


major  fraction  of  fault  planes  are  in  three  ranges:  30-60° 

(Figures  9 and  11),  uJ-bO°  (Figures  10  and  12),  and  over  30° 
(Figure  13).  In  these  cases,  waveforms  are  low  pass  filtered 
with  a corner  at  fg  to  diminish  higher  frequency  noise  such, 
as  shown  in  Figure  2. 

Half  space  notions  are  different  from  the  doubled 
whole  space  motions  for  the  dip  slip  fault  (Figures  9,  10). 

In  going  to  higher  angles  of  incidence,  the  relative  amount 
of  SV  motion  increases  and  the  whole  space  motions  agree 
less  with  the  half  space  notions.  The  static  offsets, 
determined  by  the  last  value  each  component  attains,  ;en- 


erally  differ  between  the  models.  A dislocation  in  an 

infinite  space  to  model  the  notions  in  the  half  space  uouV1 
probably  not  have  rupture  parameters  similar  to  those  used 
to  calculate  the  half  space  notions. 

For  the  strike  slip  faults  (Figures  11-13),  the 
horizontal  displacements  (including  the  static  offset)  de- 
rived from  the  whole  space  model  resemble  fairly  we1!  the 
half  space  notions,  with  the  exception  of  the  clear  "ay- 

leigh  wave  on  the  u comeoner.t  in  ™rur^  13.  At 

1 

degrees,  which  contains  loss  ?-SV  notions  than  3 f-CH  de- 
grees, the  doubled  ..hole  space  motions  approximate  t!.e 
half  space  motions  better  than  at  3 O-uO  degrees.  Above 
80  degrees,  the  Rayleigh  wave  causes  the  .agreement  to 
worsen.  The  vertical  (u^)  components  are  similar  in  "A  •- 
ure  11,  but  in  Figures  12  and  13  they  are  not.  r.”he  hori- 
zontal static  offsets  from  the  whole  space  models  match 
the  static  offsets  of  the  half  space  models.  This  ir  not 
expected  considering  the  offsets  for  point  sources  shown 
in  Figure  8,  A dislocation  model  in  an  infinite  space 
would  have  similar  rupture  parameters  to  those  used  to 
calculate  the  half  space  motions,  providing  the  verticsl 
component  was  ignored  where  appropriate. 

The  and  u^  components  in  Figure  13  show  a strong 
Rayleigh  wave,  clearly  identified  by  comparison  with  the 
whole  space  model.  Pekeris  and  Lifson  (1957)  showed  that 
for  distances  r and  source  depths  z,  the  Rayleigh  wave  is 
emerging  for  r/z  = 5 and  clearly  seen  for  r/z  = 10.  These 
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ratios  correspond  to  angles  of  incidence  of  about  79  and 
6*1  degrees  respectively.  Thus,  Figures  9-12  do  not  show 
a Rayleigh  wave. 

The  Parkfield  and  San  Fernando  earthquakes  have  been 
studied  extensively  using  high  quality,  close  distance 
accelerograms  to  derive  source  dislocation  models.  Table 
I presents  the  fraction  of  the  faults  for  these  two  cases 
which  give  direct  ray  angles  to  the  accelerograph  In  the 
ranges  9-30,  30-60,  C0-80  and  ?0-?9  degrees . "hose  per- 
centages are  estimated  for  a fault  in  a homogeneous  hnl° 
space  using  the  geometry  of  Trifunac  (197*0  for  the  fan 
Fernando  earthquake  and  of  Anderson  (197*0  for  the  r’arkfic?  ’ 
earthquake.  Even  at  the  closest  stations  to  the  fault,  at 
most  about  PO1"  o±  the  fault  is  in  the  0-30  degree  ranr*e. 

The  remaining  809  of  the  fault,  including  the  epicenters, 
is  in  the  range  where  the  free  surface  may  significantly 
distort  the  waveform  derived  from  the  whole  snace  an^roxl - 
mation.  The  San  Fernando  earthquake  had  a thrusting  mech- 
anism while  the  Parkfield  earthquake  had  a strike  slip 
mechanisr,  and  thus  S-waves  from  the  S an  Fernando  earthquake 
at  nearby  stations  would  in  general  have  a greater  component 
of  SV  type  notion.  Therefore,  when  the  infinite  space 
method  is  used,  rupture  parameters  derived  for  the  San 
Fernando  earthquake  are  morn  likely  to  be  incorrect  than 
rupture  parameters  derived  '’or  the  Parkfield  earthquake. 

Table  1 indicates  that  even  at  the  best  placed  Instru- 
ments, angles  of  inciden.ee  are  likely  to  be  high  for  meet 
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of  the  fault  plane.  This  Is  especially  Indicated  by  the 

ancles  of  Incidence  of  the  Pacoima  Dam  accelerograph.  In 

this  case,  the  Instrument  was  located  directly  above  the 

fault,  and  yet  SV-type  waves  from  only  about  20%  of  the 

fault  plane  would  be  relatively  undistorted  by  the  free 
surface . 


CONCLUSIONS 


A practical  way  of  computing  displacements  on  the  sur 
face  of  a half  space  is  by  a two  step  process.  The  first 
step  is  calculating  and  storing  the  Green's  function  for 
a given  station  for  a grid  of  points  on  the  fault;  and  the 
second  step  is  convolving  with  the  source  time  function. 
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In  studying  a particular  accelerogram  record, 
is  far  more  economical  than  it  would  be  to  recompute  the 
Green's  functions  for  each  trial  source  function.  This 
two  step  method  also  gives  a clear  idea  of  what  frequencie 
are  significant  in  the  calculated  displacement  record. 

The  dynamic  displacement  of  the  surface  of  a half  spa< 
caused  by  an  extended  dislocation  source  generally  differs 
from  two  times  the  displacement  which  would  result  from  the 
same  source  in  an  infinite  space.  There  are  two  exceptions 
The  first  exception  is  a fault  located  such  that  angle 
of  incidence  at  the  station  are  less  than  30  degrees  from 


vertical.  For  a large  tult,  however,  even  at  the  best 
Placed  instruments,  the  angles  of  incidence  will  be  this 
small  from  only  a small  fraction  of  the  fault.  This  case 
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Is  therefore  most  useful  for  studying  small  earthquakes. 

The  static  offset  estimated  from  the  infinite  space  method 
may  be  wrong  in  this  case,  but  this  is  not  of  practical 
importance  as  most  data  cannot  resolve  such  offsets. 

The  second  exception  occurs  when  the  motion  at  the 

source  is  predominantly  strike  slip.  Then  the  horizontal 
components  may  be  modeled  by  the  infinite  space  motion  when 

the  Rayleigh  wave  is  not  important  (eg.  for  angles  of  inci- 
dence less  than  80°).  The  vertical  component  in  this  case 
may  also  be  used  for  angles  of  incidence  less  than  about  60 
degrees.  This  case  will  be  most  useful  for  studying  strike- 
slip  earthquakes. 

In  many  cases,  a small  earthquake  recorded  at  a large 
angle  of  incidence  can  be  modeled  using  the  infinite  space 
method.  When  the  observed  S-wave  is  separated  from  other 
phases,  and  near  field  terms  are  small,  it  can  be  separated 
into  SV-  and  SH-  components.  Then  the  phase  shift  and 
amplitude  correction  can  be  applied  to  each  component 
separately  to  obtain  the  incident  waveforms,  for  comparison 
with  infinite  space  models. 
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APPENDIX  1 


Derivation  of  Amplification  Coefficients 
for  Body  Waves  Indicent  on  a Free  Surface 

The  ratio  of  the  amplitude  of  each  component  of  the 
free  surface  motion  to  the  amplitude  of  the  same  component 
of  an  incident  body  wave  is  derived  here.  Knopoff  et  al. 
(1957),  who  derive  the  ratio  of  each  component  to  the  total 
amplitude  of  the  incident  wave,  have  some  misprints  which 
make  this  derivation  necessary.  In  eq  6_,  a quantity  they 
refer  to  as  tann  O'  is  always  greater  than  1;  in  eq  7_,  th? 
phase  is  incorrect  in  the  center  of  the  three  equalities 
because  tan  2<J>  chances  sign  at  ^5°;  and  in  eq  9 (sub- 
stituting "tanh  0'"  as  defined),  the  phase  does  not  agree 
with  Figure  3.  The  figures  in  Knopoff  et  al.  are  correct. 
The  closely  related  equations  for  reflection  and  trans- 
mission coefficients  at  a boundary  between  two  layers  have 
a long  history  of  published  errors  (Hales  and  Roberts, 
197*0. 

We  first  discuss  the  problem  of  a plane  P-wave  inci- 
dent with  angle  3 and  a plane  SV-wave  incident  at  angle  j ; 
both  incident  from  z < 0 on  a free  surface  z = 0;  and 
both  waves  having  horizontal  slowness  s.  Later,  we  will 
specialize  to  the  case  where  only  one  of  these  waves  is 
incident . 

The  total  wave  field  is  described  by  displacement 


u = crad  4>  + curl  (o,  f,  0}  = /3$  - 34',  0,  3<{,  + 3«A 

\Dx  3z  3z"  37/ 

where 
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Here  the  first  tern  on  the  right  represents  the  incoming 
wave  in  each  case,  and  the  second  terra  represents  tie  out- 
going reflected  wave.  In  these  equations  f(t)  denotes  the 

time  dependence  of  the  potentials,  s = sin  i = sin  .1  is 

a G 

the  horizontal  slowness,  a is  the  P-wave  velocity,  and  3 
is  the  S-wave  velocity.  The  terras  P and  Sj,  are  constants, 
and  P^  and  Sr  are  unknowns  to  be  derived  from  the  free 
surface  boundary  conditions  of  zero  traction  on  z = 0. 

These  boundary  conditions  yield  the  equations: 

(pi  - Pr)  sin  21  - (S1  + S ) cos  2 j = 0 
ot2  02 

(Pi  + PrJ  C0S  2J  + (Si  “ sr)  sin  = °» 


which  have  the  solution 


p = p 
r i 


sin  2,1  - cos2  2 
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-2  sin  2J  cos  2. 


S = P 
r 


1 ^ sin  21  cos  2,1^  + Si  1 ^sln  21^ sin  2,1  - cos2  2.1^ 


where  D = 2s  cos  1 sin  2j  +_l  cos2  2j 
a g2 

The  desired  ratios  are  the  amplitude  of  a component  of  the 
motion  of  the  free  surface  (z  = 0)  due  to  either  an  inci- 
dent P-wave  (S±  = 0)  or  SV-wave  (P±  = 0)  divided  by  the 
component  of  motion  which  would  take  place  if  the 
surface  were  absent.  For  example. 


RXP  = 


ux(?s)  {-Pis  - Prs  - Sr  cos  ,J> 


up(») 

X 


-pi5 


Here  up  refers  to  the  horizontal  (x)  component  of  P-wave 
motion,  and  the  (h)  or  (“)  refer  to  the  case  of  a half 
space  and  of  an  infinite  space  respectively. 

For  incident  P-waves,  the  ratios  are: 


RXP  = 2 sin  21  cot 
a2D 


( A-l ) 


RZP  = 2 cos  2.1 
B2  D 


(A -2) 


3L* 


A-l  I 


For  S-waves  incident  at  angles  less  than  the  critical  angle 
Jc  (sin  Jc  - £),  the  ratios  are: 


RXS  = 2 cos  2,1  (A-3) 

e2  D 

RZS  *=  4s  cos  i cot  j ( A— 14 ) 


For  S-waves  incident  at  angles  greter  than  critical,  a 

phase  shift  is  introduced.  At  this  point,  it  is  necessary 

to  consider  a specific  frequency  component:  f(t)  = e^u^. 

Then  for  J > J the  quantity  cos  1 is  replaced  by  an 

a 

imaginary  quantity  — ib  sgn(u),  where  the  sign  is  choosen 
so  that  the  potential  will  decrease  with  increasing 
distance  from  the  free  surface.  Here  sgn  (w'*  = 1 if 

to  > 0 and  -1  if  u < 0,  and  b = Vs^lr7.  Then  the  ratios 

for  j > j are: 
c 


RXS 


2 cos  2 j 

e2  d' 


gi(p  + 0) 


(A-5) 


RZS 


s b cot  j 
D1 


ei(p  - £ sgn  (w)) 
2 


where 


(A-6 ) 


+ (2  s b sin 


< 


D' 


and 


L. 


A-5 


tan  p = 2g2  sb  sir,  2j  snn  (w0) 
cos2  2j 

6=1°  Jc « J 

(tt  tt/H  < J < n/2 


* *■ 
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ESTIMATED  FRACTION  (NEAREST  5 PERCENT)  OF  FAULT  AREA 
WITHIN  EACH  RANGE  OF  ANGLE  OF  INCIDENCE 


r 


FIGURE  CAPTIONS 

Figure  1:  Coordinate  system  and  faulting  parameters. 

The  plane  x = 0 is  the  free  surface.  The  f-'Milt 
strikes  in  the  x^  direction,  and  the  dip  6 Js 
measured  from  the  horizontal  plane.  The  point 

(Xj»  x2»  x 3 ) = (0»  0,  d)  is  the  geometrical 
center  of  the  fault,  and  the  origin  of  the 

coordinates  (not  shown)  which  describe 
a point  on  the  fault  plane.  The  direction  of 
slip  on  the  fault  is  given  by  the  unit  vector 
S»  wlth  components  (in  the  £ direction)  and 
S3  (in  the  C3  direction)  giving  respectively 
the  relative  amplitudes  of  strike  slip  and  dip 
slip  motion. 

Figure  2:  A model  evaluated  using  three  values  of  f , 

s 

with  the  Green's  function  in  step  1 replaced 

by  the  non-physical  function  h = 1 H(t-R/a) 

R 

(see  text).  The  table  gives  the  values  of  f 

s 

(in  Hz.),  and  the  corresponding  number  of 
time  points,  grid  points,  and  total  points  at 
which  the  function  h was  evaluated.  For  each 
calculation,  the  line  beneath  the  time  axis 
has  a length  Tg  = f“*.  For  this  calculation 
(x  , x ^ ) - (0.866  km,  0.5  km),  d = 0.5  km, 

6 s 90°.  The  rupture  model  is 

31  < 


Figure  3s 


Figure  *4: 


Figure  -3: 


Figure  6: 


Figure  7: 


a propagating  rann  with  L = 1.0  km,  = 1.0 
km,  v = 3.0  kn/sec,  t = 0.?  ace. 

The  effect  of  introducing  a 1 0 Jo  random  error 
to  each  of  the  180  evaluations  of  h (in  space 
and  time)  in  the  second  case  of  Figure  2. 
Theoretical  ratios  (solid  lines)  of  the  ampli- 
tude of  vertical  and  horizontal  components  of 
P-waves  on  a free  surface  to  the  amplitudes 
of  the  same  component  of  the  incoming  wave  (see 
appendix  equations  A-l  and  A-2)  . The  angle  of 
incidence  is  measured  from  the  normal  to  the 
surface.  The  data  points  are  these  ratios  de- 
rived from  computations  with  a point  source 
Green's  function. 

fame  as  Figure  C for  i"  •waves  (Appc rr'i::  equa- 


tions A3-AO.  The  computed  waveforms  nor  3f° 

to  80°  are  shown  in  Figure  6. 

Vertical  (u^)  and  horizontal  (u2)  components  of 
SV-wave forms  for  a half  space  (solid)  and  an 
Infinite  space  (dotted)  from  a point  source. 

For  the  source,  d = 10.0  km,  (S ^ , S ) ■ (0,  .3), 
and  6 = 90°.  The  stations  are  at  (x  , x ) = 

l 2 

(0  km,  10  tan  (angle  of  incidence)  km).  The 
source  has  a rise  time  =0.1  sec. 

Similar  to  Figure  6,  but  the  free  surface 
amplitude  and  phase  corrections  for  piano  SV- 
waves  (eqs.  A-5  and  A-6)  have  been  applied  to 

the  infinite  space  waveforms.  In  some  cases. 


the  infinite  space  waveform  has  also  been 
arbitrarily  shifted  up  or  down.  This  is  justi- 
fied because  near  field  components  are  not 
amplified  the  same  as  the  far  field  SV-waves 
compared  here.  Where  shifted,  the  dotted  lines 
which  are  disconnected  from  the  S-wave  (as  at 
80°)  show  the  zero  for  the  whole  space  waveform. 
Figure  8:  Ratio  of  static  offsets  in  a half  space  to  those 

in  a whole  space  for  a point  source  with  depth 

d = 10  km,  and  at  stations  (x  , x ) = (0  km, 

12  1 

10  tan  (angle  of  incidence)  km).  The  u cornno— 

i 

nent  is  derived  for  the  source  (£  , Z ) & 

1 3 

(5,  0);  the  u^  and  u^  components  are  derived 
for  the  source  (S  S ) = (0,  S).  "he  lines 
only  connect  points  derived  from  the  models. 
Figure  9:  Comparison  of  half  space  (solid  line)  and 

doubled  whole  space  (dotted  line)  displacements 

for  the  dip  slip  fault  with  d = 3.8  km,  (x  , x ) 

12 

= (5.0  km,  1.5  km),  <5  = 90.0°,  (3  S ) = (0,  3). 

Rupture  parameters  are  L =5.0  km,  L =3.3  km, 

1 3 

v = 3.0  km/sec,  and  t = 1.0  sec.  The  vertical 

scale  is  in  units  of  D and  the  time  scale  is 

0 

in  seconds. 

Figure  10:  Comparison  of  half  space  (solid  line)  and 

doubled  whole  space  (dotted  line)  displacements 

for  the  dip  slip  fault  with  d = 1.1  km,  (x  , x ) 

1 2 

c (5.0  km,  1.5  km),  6 = 90.0°,  (S^  S ) = (0,  5). 

Rupture  parameters  are  L = 5.0  km,  L = 1.2  km 

1 3 * 


The  vertical 


v = 3.0  km/sec,  t = 1.0  sec. 

scale  is  in  units  of  D ami  the  tine  scale  is 

o 

in  seconds. 

Figure  11:  Same  as  Figure  9,  but  here  for  a strike  slip 

fault  with  (S  , S ) = (s,  0). 

Figure  12  Same  as  Figure  10,  but  here  for  a strike  slip 
fault  with  (Sj , S ) = (S,  0) . 

Figure  13:  Half  space  (solid  line)  and  doubled  whole  space 

(dotted  line)  displacements  for  the  strike  slip 
fault  with  d = 0.5  km,  (x^  x^  = (10.0  km, 

5.**5  km),  6 = 90°,  (Sjf  SJ  = (S,  0).  Propa- 
gating ramp  model  parameters  are  L =10.0  km, 
^3  “ 0.5  km,  v = 3.0  km/sec,  t = 1.0  sec.  The 

vertical  scale  is  units  of  n and  t;  •>  tJ-o 

0 

scale  is  in  seconds. 
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8 


